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Abstract 

We consider a two-dimensional integrable Hamiltonian system with a vector and scalar 
potential in quantum mechanics. Contrary to the case of a pure scalar potential, the 
existence of a second order integral of motion does not guarantee the separation of variables 
in the Schrodinger equation. We introduce the concept of "quasiseparation of variables" 
and show that in many cases it allows us to reduce the calculation of the energy spectrum 
and wave functions to linear algebra. 
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1 Introduction 



A systematic search for integrable classical Hamiltonian systems in magnetic fields was 
started quite some time ago [HI EH HTj . A Hamiltonian containing scalar and vector 
potentials was introduced in a two-dimensional Euclidean space and additional integrals 
of motion were constructed as linear or quadratic polynomials in the momenta. The same 
problem in quantum mechanics was investigated quite recently j2j- Some new features 
have emerged in this study of vector potentials that distinguish this case from that of 
purely scalar ones. 

1. The existence of second order integrals of motion does not imply the separation of 
variables in the Hamilton-Jacobi, or Schrodinger equation. 

2. Hamiltonian systems with second order integrals of motion in classical and in quan- 
tum mechanics do not necessarily coincide 0. 

3. In quantum mechanics an additional problem arises: it is necessary to choose a 
convenient gauge in which to write and solve the Schrodinger equation. 

In the case of Hamiltonians with purely scalar potentials, quantum and classical in- 
tegrable systems may also differ, but only if the integrals of motion are third, or higher 
order in the momenta [El EH E3j ■ Moreover, third and higher order integrals of motion 
are not related to the separation of variables, at least in configuration space, in the case 
of purely scalar potentials either [9l 1 1 j . 

Variable separation in the presence of magnetic fields has also been investigated |22 
123] and turns out to be quite rare. Separable systems in this case only constitute a subset 
of the quadratically integrable ones. 

Natural questions that arise, specially in the context of quantum mechanics, are the 
following. What does one do with integrals of motion that do not lead to the separation 
of variables? How do they help to integrate the Schrodinger equation and find the energy 
spectrum? 

The purpose of this article is to provide at least a partial answer to these questions. 

The problem is formulated mathematically in Section 2. Section 3 is devoted to systems 
with one first order operator X, commuting with the Hamiltonian H. In Section 4 the 
commuting operator X is assumed to be a second order operator of a specific ("cartesian") 
type. We introduce the concept of quasiseparation of variables and identify all cases when 
quasiseparation occurs. The Schrodinger equation is solved for a specific scalar and vector 
potential in Section 5. Section 6 is devoted to conclusions. 

2 General setting 

The quantum Hamiltonian that we are considering is 

H = -y {dl + d 2 y ) - ^(Ad x + d x A + Bd y + d y B) + V, (2.1) 
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where A, B and V are functions of the coordinates x,y. The quantities of direct physical 
importance are the magnetic field Q and effective potential W: 

n = A y -B x , W = V-^(A 2 + B 2 ). (2.2) 

These quantities are gauge invariant, i.e. unchanged by the transformation 

V^V = y + (A,V0) + i(V0) 2 , 

A^A = A + V</>, A = (A,B), (2.3) 

where ^(x,y) is an arbitrary smooth function. 
The classical equations of motion are 

x = Uy-W x , y = -ttx -W y , (2.4) 

so all we need to know is Q and W (the dots indicate time derivatives). In quantum 

mechanics we need the Hamiltonian (|2.ip . Hence, given f2 and W, we must still choose a 
gauge, i.e. the function (j> in eq. (|2.3|) . and calculate A, B and V. 

The quantities Q and W are obtained from the commutativity condition 

[X,H]=0, (2.5) 

where X is the integral of motion, i.e. either a first, or a second order linear differential 
operator. The integral X is obtained from the same condition. 

The two equations to solve simultaneously, once a gauge is chosen, are 

H?P = EiP, Xrp = Xip, (2.6) 

and the vector potential (i.e the gauge) should be chosen so as to simplify this pair. 

3 First order integrals 

A first order integral of motion will have the form j2] 

X = a(L 3 + yA- xB) + 0(P X + A) + 7 (P 2 + B) + m, (3.1) 

where L3, Pi and P2 are the angular and linear momentum operators, i.e. 

L 3 = -ih(yd x - xd y ), Pi = -ihd x , P 2 = -ihd y , (3.2) 

a, P and 7 are constants and the functions m(x,y), £1 and W satisfy 

(ax - 7)0 + m x = 0, (ay + 0)tl + m y = 0, (3.3) 
(ay + (3)W X + (-ax + j)W y = 0. 
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Two inequivalent possibilities occur. For a we can put a=l and translate (3 and 7 
into /3=7=0. For a=0, /3 2 +7 2 7^0 we can rotate (3 into /3=0 and normalize 7 to 7=1. Let 
us now look at the two cases separately and solve the Schrodinger equation in each case. 

a) a=0, /3=0, 7=1 

We have [2j 

m = m(x), W = W(x), O = U(x) = rh(x). (3.4) 
The integral of motion (|3.1|) in this case reduces to 

X = P 2 + m{x)+B{x,y). (3.5) 
A convenient choice of gauge is 

B = -m(x), A = 0, (3.6) 

and from (12.211 we obtain 

Q = rh(x), V = W(x) + ^m(x) 2 . (3.7) 

The system (|2,6|) reduces to 

h 2 1 
{" y(S£ + Sj) + ihm(x)d y + W(x) + -m(x) 2 }^ = 

= (3.8) 
With this choice of gauge we have the separation of variables, i.e. 

Hx,y) = f Ek (x)e-f y , (3.9) 

h 2 d 2 fk 2 1 

~ Td^ + [ T + fcm + w + 2™* ~ E]f = °' (3 ' 10) 

Equation l|3.1()|) is exactly solvable for many choices of m(x) (i.e. for many choices of 
the magnetic field Q) and the scalar potential V(x). For instance, if we have 

m = u 2 x 2 , n = 2co 2 x, V(x) = 0, (3.11) 

eq. (|3.10|1 has the solution 

f = e-^-H n (rx), T = J^(2k) 1 '\ E = + (2n + l)^-V2k, (3.12) 
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where H n (£) is a Hermite polynomial. In this case f(x) is regular and square integrable, 
however the solution l|3.9p involves a plane wave in y. 
b) a=l, /3=7=0 

In this case the determining equations (3.3) imply [2] 

m = m(r), W = W(r), fi(r) = , (3.13) 

r 

where we are using polar coordinates 

x = rcosQ, y = rsinQ. (3-14) 
The operator (|3 . l|i in this case is 

X = L 3 + yA- xB + m, (3.15) 
and a good choice of gauge is given by 

yA- xB + m = 0, (3.16) 

leading to 

A( r , 6) = - r ^-sin@, B(r, 6) = ^-cosQ, (3.17) 
r r 

2 

fry 

V(r,Q) = W(r) + - 1 . (3.18) 
The solution of the system l|2.6jl can be written in separated forms as 

^(r,e) = e- me ^R E>M (r), (3.19) 



with R(r) satisfying 

h 2 d 2 R rrTT , , 1 r h 2 M 2 m(r) 2 _ h 2 „ ^„ 

~ y^ 2 + {w{r) + V 2{ ^r~ + 2^ " hm{r)M ~ ~8 ]}R = K (3 0) 

Again, this equation is exactly solvable in special cases, e.g. 

, s ol a 2 , s 

m(r) = -, W = -— 4 . (3.21) 

The conclusion from this section is that first order integrability in a magnetic field 
leads to a separation of variables either in cartesian, or in polar coordinates. To make this 
happen, a proper choice of gauge is crucial. Indeed, in a previous article [2] a different 
choice of gauge was made, leading to R-separation, rather than ordinary separation. The 
separated ordinary differential equations (|3,l()p and l|3.2()|l both have the form of one- 
dimensional Schrodinger equations. The magnetic field and the effective potential combine 
together into an x-dependent, or respectively r-dependent one-dimensional "potential". 
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4 Second order Cartesian integrability and the qua- 
siseparation of variables 

Let us now consider a Hamiltonian of the form 112.111 . admitting one second order integral 
of motion X. It has been shown that this operator X (or function in classical mechanics) 
can be transformed into one of four standard forms [2,3]. They were called Cartesian, 
polar, parabolic and elliptic, because in the absence of a magnetic field, their existence 
leads to the separation of variables in Cartesian, polar, parabolic or elliptic coordinates, 
respectively [HlEO]- The direct relation to the separation of variables in the Schrodinger, 
or Hamilton-Jacobi equation does not hold in the presence of a magnetic field, but we keep 
the terminology. In this article we restrict ourselves to the Cartesian case. 
The "Cartesian" integral of motion has the form (21 El 

X = -yd x 2 - ih[(A + h)d x + k 2 d y ] - j(A x + k lx + k 2y ) + X -A 2 (4.1) 

+m + k\A + k 2 B. 

All functions involved in the Hamiltonian (|2,1|) and the integral Il4.ll! can be expressed 
in terms of two functions of one variable each, / = f{x) and g = g(y), satisfying 

f = af 2 + (3f + 7 , g" = -ag 2 + 5g + £, (4.2) 

where a, (3, 7, 5 and £ are real constants. We shall also use the first integrals of eq. (|4.2|l . 
namely 

f 2 = Ps(f), g' 2 = Qs(g) (4.3) 

W) = l a f + + 2 7/ + <*U Q-M = -\ag* + Sg 2 + 2£g + a 2 

where a\ and a 2 are further real constants. The dots and primes are x and y derivatives, 
respectively. 

Eq. (|4.2|) and l|4.3jl can be solved in terms of elliptic functions, or their degeneracies, 
if the polynomials P${f), or Qz{g) have multiple roots. In terms of f(x) and g(y) we have 

n = f(x)+g"(y), 
W = -^f-gf-^(f-g) 2 + (Z- 1 + v)(f-g), 

ki = -g'(y), k 2 = -f(x), (4.4) 
m = -f (g 3 + 2f - 3gf 2 ) + 0{fg - f 2 ) - & -{f 2 - g 2 ) 
->y(2f-g)+nf + Zg. 
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where fi, figuring in the effective potential W and in m is an additional constant. The re- 
sults 114.211 and (4.4) were obtained in classical mechanics [3], but the classical and quantum 
results coincide in the case of Cartesian integrability [2]. 

The vector potential (A,B) is yet to be chosen, but must satisfy 

A y -B x = VL = f{x)+g"{y). (4.5) 

We shall in this section assume f 2 + g" 2 ^ 0. 

Thus, we have a Hamiltonian H and first integral X satisfying (12.111 and (|4,1|) . respec- 
tively, expressed in terms of quantities satisfying eq. I|4.4p . In general, variables do not 
separate in eq. I|2.6|) in any system of coordinates. 

Let us introduce the concept of "quasiseparation of variables". We take a linear com- 
bination of the two equations (|2.6|l . namely 

{(H-E) + <f>(x,y){X-\)}V = (4.6) 

where 4>(x, y) is a function to be determined. We wish to choose the function (fi(x, y) in 
such a manner that eq. (|4.6|) allows separation in Cartesian coordinates. The solutions of 
eq. (|4,6|) will then have the separated form 

^e\{x, y) = v EX (x)w EX (y). (4.7) 

The solution of the Schrodinger equation will be a superposition of separated solutions: 

■® E (x,y) = J A EX VEx(x)wEx(y)dfi(X) (4.8) 

where /i(A) is some measure to be chosen and A E \ is independent of x and y. For bound 
states the integral in eq. I|4.8p will reduce to a sum. 

On a more formal level we introduce the following definition. 

Definition 1. The commuting pair of operators {H, X} allows the quasiseparation of 

variables in the system 12. 6)) if there exists a function <fi(x, y) such that eq. \4.6\ j allows 
the separation of variables in the sense of eq. j4- T\) . 

In this article we are considering the case when (x,y) are Cartesian coordinates in a 
plane, but the concept of quasiseparation is easily generalized to other coordinates and 
other spaces. 

Substituting (|4,7|) into lj4.6p we obtain the equation 

^ vw + <ftlVW + (f>2VW + (f>3VW + <fi/±VW = (4.9) 



with 

h 2 



, > x l + 0), ^2 = -ihA{l + 4>)+ih(j)g , fo = -iK{B-4>f), (4.10) 
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'-{{1 + 4>)A X + B y )-E-X(j). 



fa = W + ^A\l + 0) + l -B 2 + 0(m - g'A - fB) (4.11) 

ih . 

The necessary and sufficient condition for variables to separate in eq. I)4.9p is that we 
have 

01 V 1 (x) cp 2 V 2 {x) 



04 V{x) + W{y) 04 V(x) + W(y) 
-h 2 W x {y) 03 W 2 (y) 



(4.12) 



204 V(x) + W(y) 04 V(x) + W{y) 

where V%, Wi, V and W are some functions. 

Let us consider the choice 4>(x,y) = —1 separately. 



4.1 (f)(x,y) = -l 

With this particular choice, eq. 1)4.9)1 simplifies, V^x) reduces to a constant, and we find 
that 

B(x,y)=-f + ^^-g, and W x {y) = (4.13) 
V 2 2i g 

Condition 1)4.5)1 becomes 

A(x,y) = g+r(x), (4.14) 



where r(x) is arbitrary, and eqs. (|4. 12(1 reduce to 

9 V 2 



04 ih(V(x) + W(y)) 

Thus 



(4.15) 



<)2 -&)=0 (4.16) 



is a necessary and sufficient condition for quasiseparation. A straightforward calculation 
shows that this condition only depends on g(y) and is equivalent to 

fin I fit 

9 2 -9 9 = 0. (4.17) 

The general solution of eq. ()4.17)l is g(y) = C\e c ' 2V + C3 and a two parameter class 
of particular solutions is g{y) = g\y + go. It follows that c/>(x,y) = —1 only allows 

quasiseparation in the cases a = 5 = £ = and a = 0, 5 > (i.e. <?(y) = g\e^ y + 
g2 e'V5y + |) with gl52 = o. 
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4.2 (j)(x,y)^-l 

Eliminating (f>4 from the equations (|4. 12|l and using relations (|4- 1 0|1 and H4.11|) we obtain 

A(x, y) = r{x) + ^ g (y), B(x, y) = V (y) + 0(s, y)f(x), (4.18) 
1 + (f)(x,y) 



e 2 ^c 2 + 2kg( y )^ f(x) 2 ' 



T /f \ / ci - 2kf(x) c 2 + 2kg(y) 

Fl( )=ei V /(-) 2 ' l(y) = e2 V^)^ ( ] 

V 2 (x) = jr^V^x), W 2 (y) = jv(y)W 1 (y). 

The functions r(x) and r/(y) are arbitrary and can be modified by gauge transformations 
(for instance they can be set equal to zero). The entries c\, c 2 , k, e\ and e 2 are constants 
with ef = ef = 1. 

The conditions 1|4.18|) on A, B and 4> are necessary for separation of variables in eq. 
1|4.9|) . There is a further necessary condition that together with (|4. 18|l is sufficient, namely 

W 1 (y)fa(x,y)=X(x) + Y{y), (4.20) 

where X(x) and Y(y) are arbitrary functions. 

In other words we must determine the conditions on the functions f{x), g{y), t{x), 
7](y) and the constants k, ci, c 2 , e\, e 2 in such a manner that eq. (|4.20p is satisfied. All 
of the above quantities are real and the square roots \Jc\— 2k f and \Jc 2 + 2kg must be 
simultaneously real, or simultaneously imaginary. 

To proceed further we use l|4.11jl . I|4.18|) and (|4.19jl to obtain 

Wt(y)Mx, V) = R(x) + S(y) + T(x, y) (4.21) 



R(x) = 

c 2 f 2 (x) + ( Cl - 2kf(x))(2\ - 2a x + a 2 - 2/x/(x) + 5f 2 (x) - t 2 (x) + ihf{x)) 



2e ly / Cl -2kf(x)Jf(x) 



(4.22) 



S(y) = 

cig' 2 (y) + (c 2 + 2kg(y))(2e - 2\ + a x - 2a 2 + 2 m (y) + Pg 2 (y) - rfjy) + ihrj (y)) 

2e 2 y/c 2 + 2kg(y)y/g'(y) 2 

(4.23) 
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and 



T {x iy )= F ^a(y) + f( x ) G i(y) 



v/ci - 2kf(x) V c 2 + 2%(y) 



^(^W^^ (4.24) 



F = ^[-^ 2 (^) + ( / ^ fc/( x))/(x)] i 



F 2 (x) 
Gi(y) = 
G 2 (y) ^ 



In (|4.24|) we have 

. e 2 [kf 2 (x) + ( Cl -2kf(x))f(x)] 

2( Cl -2kf(x))f(x) 
e2[kg' 2 (y) + (c 2 + 2kg(y))g"(y)] 

e 1 [kg' 2 (y)-(c 2 + 2kg(y))g"(y)] 
2{c 2 + 2kg{y))g'{y) 

The separability condition l|4.2()|) is equivalent to the condition 

d 2 T 



dxdy 



0. (4.26) 



We can consider the real and imaginary parts of eq. I|4.26j) separately. We obtain two 
conditions: 



1 / Ft v 1 f °i 



f\y/c 1 -2hfJ g'\^c 2 + 2kg 



Ni (4.27) 



! C1 - 2k f _ kC > c 2 + 2kg 
P 1 9 



kF^r-^± = kG'^r^jf* = N 2 , (4.28) 



where Ni and N 2 are constants. 

More explicitly, eq. (|4.27|l and l|4.28|) can be rewritten as 

- gl [fc 2 / 4 + 2fc( Cl - 2kf)f 2 f + ( Cl - 2fc/) 2 (/ 2 - //)] 

(d-2fe/)3/2(/2)3/2 

6 2 [fc 2 g 4 - 2fc(c 2 + 2%)gV + (eg + 2%) 2 (/ 2 - g'g'")} 
{c 2 +2kgf/ 2 {g' 2 f/ 2 



= Ni (4.29a) 
Nx (4.29b) 
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ke 2 [2k 2 f* + fc(ci - 2kf)pf - (ci - 2fc/) 2 (/ 2 - //)] = 2g 
2( Cl - 2fc/)3/2/3 2 

-ketftfg' 4 - fc(c 2 + 2fcg) g y - (c 2 + 2fcg) 2 (g" 2 - g'/')] 

2(c 2 + 2%)3/V3 = N2 (429d) 

We see that the functions r(x) and r](y) do not figure in ( |4.29j) . and hence have no influence 
on the separation of variables. The functions f(x) and g{y) depend on the constants a, /3, 
7, 5 and £ of eq. 1)4 .2|) and on a total of four further integration constants. Our aim now 
is to find all values of these constants and of k, c\ and c 2 for which eq. 1)4.29)1 are satisfied. 
Let us consider the cases k = and k ^ separately. 

I) k = 

Prom (l4~2H we have N 2 = and H4.29cl4.2ydB are satisfied identically. Eq. H4.29al4.2yhH 

simplify to 

^ = ei ^ r " ( j2 ) 3/2 = e2 ^ g (^)3/2 • ( 4 - 30 ) 

Using eq. 1)4 .2|l we obtain 

3eiATi/ = 2aVHr// 2 , 3e 2 W = 2a^^. (4.31) 

Eq. ()4.31)l are only compatible with 1)4.2)1 if we have a = N\ =0. The assumptions ci = 
or c 2 = lead to contradictions or <ft(x,y) = —1, so we are left with f 2 — /(/) = 0, 
9 2 — 9 {g )=0 an d hence, using 1)4.2)1 again we find that the only solutions for f(x) and 
g(y) are 

/(x) = / ie v^ + /2e -v^_2 

<7(l/) = 5ie^ + g 2 e-^-| (4.32) 
/5>0,5>0,/i/ 2 = gig 2 = 0, 

where f{ and gj are constants, or one of the functions f(x) or g(y) may be linear, the other 
being as in 1)4.32)1 . 

II) k ^ 

We shall run through different possible solutions of eq. 1)4.2)1 and determine which of 
them are compatible with eq. 1)4.29)1 . 

1) a = f3 = 8 = 0, 7£^0. 
We obtain 

f(x) = ]px 2 g(y) = ^y 2 (4.33) 
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Eq. igjH) and the condition (c 2 + 2kg(y))(c 1 - 2kf(x)) > imply 



Cl = c 2 = iVi = iV 2 = 0, £7 < (4.34) 

2) a = 0, /? 2 + <5 2 / 0. 

The functions f(x) and g(y) will be expressed in terms of exponentials, trigonometric 
functions, or one of them may have the form 1|4.33|) . In none of these cases can eq. (|4.29|l 
be satisfied. 

3) a ^ 0. 

In this case we are dealing with the two nonlinear equations (|4.3|) . If the cubic polyno- 
mial on the right hand side has three distinct roots, we obtain solutions in terms of elliptic 
functions. Otherwise the solutions involve elementary functions. 

Let us discuss the equation for f (x) . We can assume with no loss of generality that we 
have a > 0. Indeed, if we replace f(x) — > —f{x),a — ► —a, (3 — > /3, 7 — > —7,0"! — > o\ in 
eq. I|4.3p we get the same equation. Hence we can change the sign of a from negative to 
positive (if necessary). We write 

f(x) 2 = ~a(/(s) - h)(f(x) - f 2 )(f(x) - / 3 ), a > (4.35) 

We require f(x) to be real, hence f(x) 2 > 0. If all three roots are real we order them as 
fl < fi < /3- Otherwise we consider fx G M. and f2=p + iq,f3=p — iq,p,q£^,q> 0. 
The possibilities are: 



4 /(*) = -£-, 

2a ' 2a a(x — xq) 



a) f 1 = f 2 = f 3 = -£-, f( x ) = -£- + — (4.36) 



&) /i = / 2 </ 3 </, /(a) = /i+ . 2 /S , A v (4.37) 

Sin W(X — Xq) 



-/? + 2^ -4a 7 2 y 7 /? 2 - 4a 7 6^ 2 
/l = Ya ' U = 4 ' (/3 " /l) = — 



h<f<h = h, f(x) 



(3 + 2 V/? 2 - 4a 7 



2a 



a; 



/1 + (/3 - /1) tanh u{x - x ), 

6co 2 
a 



VP 2 - 4«7 



(/3-/l) = 



(4.38) 



d) fi<f2 = h< f, f{x) =h + 



/3 - /l 



tanh 2 uj{x — xq) ' 



(4.39) 
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with fi, / 3 and ui as in (|4.38|l . In cases b,c and d we assume /3 2 — 4«7 > 0. 



e) h<f<h<h,h<f-2 fix) = fi + (f 2 -fi)sn 2 (u(x-x ),k), (4.40) 



2 «/r j- >, //2-/l 



^ = ^(/3-/l), fc 



where sn(o;a;, fe) is a Jacobi elliptic function. 

/) A</2</3</, f{x) = h+ 2 /f~ h s hV (4-41) 
with k and u; as in 114.4 



5) /l€R, f 2>3 =P±iq, q>0, f( x ) = f l+ A 1Cn(pX ^ (4.42) 

1 + cn(px, k) 

/2a 2 2 2 2 

P= VT = — 2I — ' =(p-A) +9 

The solutions for g(y) are similar and can be obtained from those for f(x) by the 
substitutions 

f{x) -> ~g(y),a -> -a, /3 -> (5, 7 -> 01 -> cr 2 (4.43) 

The solutions for /(x) and g(y) can now be substituted into eq. (gHH) in order to 
determine whether there exist constants c±, c 2 , k and e±, e 2 for which the quantities N\ 
and N 2 are indeed constant. It turns out that if f(x) or g(y) are given by elliptic functions, 
eq. (|4,29p are never satisfied. However, if both polynomials in eq. l|4.Hp have multiple 
roots, eq. (|4.29p can always be satisfied. We give the results for f(x) and g(y) in Tables 
1 and 2, respectively. The calculations are quite cumbersome and were performed using 
Mathematica. 

Prom the tables, we see that any f(x) from Table 1 can be combined with any g(y) 
from Table 2. The fact that N 2 must be the same in both tables provides a relationship 
between e\ and e 2 . 

More specifically, we have: 

ei = e 2 for (iq, G 3 ), {F 2 , G 3 ), {F 3 , G x ), (F 3 , G 2 ), (4.44) 

(F 3 , G4), (F4, G 3 ) 

ei = -e 2 for {F 1 ,G 1 ),(F 1 ,G 2 ),(F 1 ,G 4 ),(F 2 ,G 1 ),(F 2 ,G 2 ), (4.45) 
(F 2 , G4), (F 3 , G 3 ), (F 4 , Gi), (F4, G 2 ), (F4, G4) 
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Similar tables are easily obtained for a < and we shall not present them here. We 
see that for a ^ we must have k ^ (k = would imply c\ = c% = 0), Otherwise, k 
remains arbitrary, as do a, in eq. (|4- 2f) and 1)4 .3|) . The integration constants 0"i and 
o"2 must be such that the polynomials Ps(f) and Qs{g) in eq. (|4. 3|) have multiple roots. 
The constants c\ and C2 are completely determined. 

Let us sum up our results clS ct theorem. 

Theorem 1. Separation of variables in eq. \4-b]) /or f 2 +g 2 ^ occurs if and only if we 
have one of the following 

1. a = 5 = (,=0. 

In this case we have g{y) = goy + gi and f{x) is any solution of eq. \4-.2\) with a = 0. 
We can put 4> = — 1, A = g and B = — /. 

2. a = 0, 5 > 0, 

g(y) = gi e^y + g 2 e~^y - J, 5 i<?2 = 

and /(x) is any solution of eq. \4-2\j with a = 0. 
We can pw£ = —1, A = g and B = —f. 

3. a = (3 = 5 = 0, £7 < 0, a = c 2 = 0, k + 0, (f> = y^. 
In this case we have f(x) = \^x 2 , g(y) = \^y 2 
and we can put A = j^rg and B = <pf. 

4- ct 7^ 0, k 7^ 0, fix) and g(y) are solutions of eq. \4-2\) as listed in Table 1 and Table 
2, respectively, and <\>, A and B as in eq. \4- 18§ where we can put r = n = 0. The 
values of the constants c\ and c 2 are listed in the tables and e% and e 2 are related as 
in eq. \4-44\) or \4-4$) - 

The formulas for the vector potential (A, B) can be modified by putting A — > A+t(x), 
B — > B + rj(y) without any effect on solutions. 

5 Example of solving the Schrodinger equation by 
quasiseparation of variables 

In order to show how quasiseparation of variables allows us to solve the Schrodinger equa- 
tion, let us consider the simplest case, namely (|4.33|) . We have f(x) = g(y) = 
and hence 

Q = Q = 7 + e, W = i(e - 7 + V){ix 2 - £y 2 ), 7 £ < (5.1) 
Changing the notation of the constants, we put 

W = ^{uolx 2 + uo 2 y 2 ), n = n , ui^uj 2 , (5.2) 
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In eq. 1)4.18)1 we choose t(x) = rj(y) = and obtain 



-y, B = ■ x, <p{x,y) = . (5.3) 



In terms of the constant magnetic field fto and the frequencies uj\ and u>2 the Hamil- 
tonian H and integral X reduce to 

H = -%(d 2 x + dl) - iH Q ° fayd x - uj lX dy) + (5.4) 

l fa+0J 2 ) 2 + ^Q f 2 2, 2 2^ 

2 fa+^y ^ + UJ * y } ' 

V ^ «2 * ^0^1 1 2 (Wi + W 2 ) 2 + S^o 2 , K n 

X = — -5 - ih^ ofaydx - toixdy) + -uj x — ■ yr^-x . (5.5) 

2 uf — U2 2 (LU1+UJ2) 

The linear combination (|4.6|l of H and X that allows the separation of variables is in 
this case 



{-% l-dl + dl] + ] [fa + u 2 ? + nljfa^ + Wa y )^ 

2 CJl J 2(Wl + CJ2 ) 



_ ^^ A) ^. (5 . 6) 



Putting 

ip(x,y) =v(x)w(y), (5.7) 

as in eq. 1)4 .7|) we obtain 

, 1 ^1^2 „ . ^2 , n 2i 2 . ^1 ~ ^2 , , , K „s 

- — V xx + -- ; ~^[fa+UJ2) +U \xv-\ Xv = k v, (5.8) 

2 Ml 2 [Ui + CJ2) &l 

^-2 ^2 

—?T w yy + 7T7 r — ^2 K^i + w 2) 2 + ^o]y 2 ^ - Ew = -k w, 

l Z(UJi + U>2) 

where ko is a separation constant and (E,A) are fixed constants. The solutions of eq. 1)5.8)1 
are of course well known and the regular, square integrable solutions are given in terms of 
Hermite polynomials H n as 

v(x)=e 2 H ni (Tix), w(y) = e 2 H n2 (T 2 y). (5.9) 
The constants satisfy 

' ] "° , [(o; 1 + ^2) 2 + ^] 1/4 , « = 1,2, (5.10) 
nfa + W2J 
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k LUi — (u)\ — UJ 2 )X = 

E-k = 

n i= 0,l,2,... 

Eliminating ko we have 

Eu>l + \{uJ 2 — ^>l) 

n 



2{tj)\ + u 2 ) 
2(ui + u 2 ) 



hw\uj 2 

LOl + Ul 2 

n\ + n 2 . 



(wi+£j 2 ) 2 + fi§(2ni + l) 



(cJi+w 2 ) 2 + ^(2n 2 + l), 



(5.11) 



{uji +uj 2 ) 2 + Ql(n+ 1) 



(5.12) 



The solution ip(x ; y) thus depends on two nonnegative integers n\ and n 2 , whereas the 
constants E and A depend only on n. 

Let us now fix E and A, i.e. fix n = n\ + n 2 and write the solution of the Schrodinger 
equation with Hamiltonian (|5.4|l as a superposition of solutions of eq. Q5.6p : 



e 2 



^ ^n 1 ,n-ni-ffn 1 (Tlx)i? n _ ni (r 2 y), 



(5.13) 



ni=0 



where the constants ^4 nijn _ ni are to be determined. We substitute 1)5. 13j) into the Schrodinger 
equation, use the recursion relations for the Hermite polynomials and obtain a linear ho- 
mogeneous equation for the constants A ni ^ n _ ni . This equation is best written in matrix 
form involving a tridiagonal matrix 



M\A >= |0 >, 



(5.14) 



i.e. 



/ an —S 
nS a 22 — 2S 




W A 0n \ ( ^ 



V 





25 a nn —nS 
S a n+ i n+ i J 



A 



ln-l 



A n -n 
V A n0 j 



V / 



(5.15) 



an = R(wi + (2n + 1)uj 2 ) - E, a 22 = i?(3wi + (2n - l)w 2 ) - E, 
= R((2n - l)wi + 3lo 2 ) - E, a„+i n+ i = R((2n + 1)^ + uj 2 ) - E, 

R _ fiyVi + u; 2 ) 2 + n% s _ ihn 

2(oji + ui 2 ) u)\ + a; 2 



(5.16) 



■yJu)\U) 2 . 
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The energy is obtained from an algebraic equation of order n+1, namely 

detM = 0. (5.17) 

The constant A, for E and n given, is obtained from eq. (|5.12|l . The wave function for E 
and A fixed has the form (|5.LSp with coefficients A nin _ ni obtained by solving the system 
<|5T5ll . 

Thus, once eq. I|5.6|l is solved by separation of variables, the solution of the Schrodinger 
equation reduces to linear algebra. The energy operator is reduced to block diagonal form, 
with each block of finite dimension, namely n+1. In this scheme the problem is "exactly 
solvable". 

For low values of n this can be done explicitly. For the ground state we have: 

h r, 7T7 . 



n = E = -J( Ul+uj2 )2 + nl A = - ^—J^ + ^ + nl (5.18) 

Z v + U>2) v 

rfx 2 +T$y 2 

4>oo = A 00 e 2 
The first two excited states satisfy 
h 



n = l E± = -[2^{uj x + co 2 ) 2 + Q 2 ± yjiux - u 2 ) 2 + fig], 

h 



A ± = 2t[-M V^l +"2? + fig (5.19) 

Z{U 2 — UJ l ) v 



+ u) 2 )\J {oji - uj 2 ) 2 + Qq], 
Tp± = e~ ' 2 *~ [A^Hifoy) + A^i(tu)], 



2 2 2 2 



with 



Af = " 2^[y + ^) 2 + + 3oj 2 ) ±yj{u x - + + u 2 )}A± v (5.20) 

For n=2 and n=3 we must solve a cubic and quartic equation, respectively. 



6 Conclusions 

We have shown that in the presence of a magnetic field VL ^ the existence of a first order 
integral of motion leads to the separation of variables in the Schrodinger equation, once 
we make the proper choice of coordinates and gauge. 

The existence of a second order integral does no longer imply the separation of variables. 
In many cases, identified in Theorem 1 of Section 4, we can '"quasiseparate" variables, 
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that is separate variables in one equation that is an appropriate linear combination of the 
two equations, Hip = Eifi and Xip = Xip. The solutions of the Schrodinger equation are 
then linear combinations of solutions of this equation, namely (|4.6|) . 

In Section 5 we analysed one of the separable cases in detail, namely that of a constant 
magnetic field Qq and an anisotropic harmonic oscillator effective potential W (see eq. 
Q). 

Thus quasiseparation of variables is reminescent of the Dirac equation [121 • The same 
method can be used for all cases identified in Section 4. 

The question that remains is the following: What can we do in those integrable cases 
when neither separation of variables, nor quasiseparation of variables occurs? The same 
questions arise in the case of polar integrability where the magnetic fields and effective po- 
tentials are known [21 E], but the Schrodinger equation remains to be solved. Some exam- 
ples of parabolic and elliptic integrability are known the solutions of the Schrodinger 
equation remain to be studied. 

Finally, a few words about superintegrability [U EI El EI [113 CH [T21 [T^J [TH1 EE! [IH1 HH1 
123 EI] In two dimensions, superintegrability means that two independent operators, X\ 
and X%, commuting with H (but not with each other) exist. The only known case with 
a nonzero magnetic field that exists, is that of a constant magnetic field fio and a zero 
effective scalar potential W=0. In this case there are three first order integrals and they 
generate a four- dimensional Lie algebra, isomorphic to a central extension of the Euclidean 
Lie algebra e(2) [2j. 
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Table 1 

Functions f(x) leading to quasiseparation of variables. We have 

u x = , /3 - 4a7 > 0, a > 0, / 3 - /i = , Ni = — 

4 a kt\ 



No. 


m 


Cl 


iV 2 


Comment 


Fi 


6 /3 
a(x— xq) 2 2a 


2fc/i 


/ —ka 


h = h = h< f(x) 
f - P 

J 1 2a 


F 2 


f 1 /3-/1 

SUl 2 U)\(x — X{)) 


2A;/ 3 


r k 1 —ka 
^2\j — 


h = h<h< f(x) 

f -/3-v//3 2 -4a 7 
J 1 2a 


^3 


fi + if 3 - /i)tanh 2 u x (x - x ) 




r k —ka 
~^2\J — 


fi < f{x) <h = h 

^ 3 2a 


F 4 


f 1 /3-/1 

•'l tanh 2 u>i(x— xq) 


2fc/i 


r k —ka 
^2\j — 


h < / 2 = h < fix) 
, _ -p+^Jp 2 -^ 

■> 3 2a 



Table 2 

Functions g(y) leading to quasiseparation of variables. We have 
V^ 2 + 4af 



,S 2 + Aa£ > 0,a > 0,g 3 - g x 



6u4 
a 



e 2 N 2 



No. 



g(y) 



c 2 



Comment 



6 



2a a(y-y ) 2 



-2kg 3 



* k 



-ka 
3 



g(y) <g\ = 92 
93 



93 



G 2 



n 93-91 

y3 + Q „vi2 , 



tanh 2 ui 2 {y-yo) 



-2kg 3 



-ei- 



-ka 
3 



g(y) < 91 = 92 < gs 

<5- v /<5 2 +4ag 
Pi - 2a 



5-3 - (P3 - Pi)tanh 2 u; 2 (y - y ) 



-2kg 3 



61 2 



-fcq 
3 



Pi = #2 < p(p) < P3 



Pi 



2a 



G*4 



„ _ 93-91 
^ 3 sin 2 W2(y-yo) 



-2k 9l 



c k 
" 6 l2 



-ka 
3 



p(p) <9i< 92 = 93 

_ 5+^/lfl+Aa~i 
93 - 2a 
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